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Abstract 

Geometrical aspects of quantum computing are reviewed elementarily for non- 
' experts and/or graduate students who are interested in both Geometry and Quan- 

O ■ turn Computation. 

In the first half we show how to treat Grassmann manifolds which are very 
important examples of manifolds in Mathematics and Physics. Some of their ap- 
plications to Quantum Computation and its efficiency problems are shown in the 
second half. An interesting current topic of Holonomic Quantum Computation is 
£h I also covered. 

In the Appendix some related advanced topics are discussed. 
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1 Introduction 

This is a review article based on lectures given at several universities in Japan and a talk 
at Numazu-meeting[|. The aim is to show a somewhat unconventional but fruitful path 
connecting Geometry and Quantum Computation, and the audience is graduate students 
and/or non-experts who are interested in both of the disciplines. 



The progress of Quantum Computation after the excellent work of P. Shor [2S] on prime 



factorization of integers and the work of L. Grover |L2[ on quantum data-base searching 
is very remarkable. These discoveries have given great impacts on scientists. They drove 
not only theoreticians to finding other quantum algorithms, but also experimentalists to 
building practical quantum computers. For standard introduction see for example f2lfl . 

\m or ma. 



The conventional methods of Quantum Computation are more or less algebraic. On 
the other hand we are interested in geometrical or topological methods. Geometry or 



*E-mail address: fujii@yokohama-cu.ac.jp 

1 A meeting held by Yoshinori Machida at Numazu College of Technology to discuss recent results on 
Geometry, Mathematical Physics, String Theory, Quantum Computation, etc. 



Topology are crucial to understand mathematical or physical objects from the global 
point of view. 



For general introduction of Geometry and Topology the book [18] is strongly rec- 
ommended. But in this book the volume calculations of some important manifolds like 
Grassmann ones or more generally symmetric spaces are missing. They are important in 
understanding of entanglements or entangled measures. In the first half we show in some 
detail the volume calculations of Grassmann manifolds. Here let us recall some basic 
concepts. 

A homogeneous space is defined by 

M £ G/H, 

where G is a Lie group and H its subgroup. We are particularly interested in the case 
where G is a classical group (for example, a unitary group U (n) or an orthogonal group 
0(n)). The complex Grassmann manifold Gk, n {C), which is our main concern in this 
paper, is written as 

hA } ~ U(k)xU(n-k)- 

The volume of Gk, n is expressed in terms of the well-known volume of U(n) (see the 
following sections), 

Vol(« M (C))- 



Vol ([/(£;)) x Vol (U(n — k))' 



This is the usual method to obtain the volume of homogeneous spaces. 

On the other hand, the volume is obtained by integrating the volume form of Grass- 
mann manifolds (: dv(Z, Z*)) which is expressed in terms of local coordinates (: Z) (see 
the following sections): 



Vol(G fc) „(C)) = / dv(Z,tf) 



Is it really possible (practical) to carry out the integral on the right hand side? As far 
as we know such a calculation has not been performed except for k — 1 (case of complex 
projective spaces). For k > 2 direct calculation seems to be very complicated. We would 
like to present this calculation as a challenging problem to the readers. 

Let us come back to Quantum Computation (QC briefly). 
Gauge theories are widely recognized as the basic ingredients of quantum field theories 
which enjoy remarkable progress recently, String Theory, M-Theory, F-Theory, etc. There- 
fore it is very natural to incorporate gauge theoretical ideas to QC; that is construction 
of "gauge theoretical" quantum computation and/or of "geometric" quantum com- 
putation in our terminology. The merit of geometric (or topological) method of QC may 
be the stability with respect to the influence from the environment. 

In |3]J and ]22| Zanardi and Rasetti proposed an attractive idea • • • Holonomic 



Quantum Computation • • • using the non-abelian Berry phase (quantum holonomy in 
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the mathematical terminology). We introduce this concept in the final section. See also 
16] and |24| for another interesting geometric model. 

Quantum Computation comprises of many subjects. To give a comprehensive overview 
is beyond the scope of this article, so we focus our attention on the construction and the 
efficiency of unitary operations, and give geometric interpretation to them. Here let us 
make a brief review. 

For n = 2* (t G N) we set a unitary operation 

U f : (C 2 )®* — > (C 2 )®' ; U f (\a)) = (-l) /(a) |a) 
where / is a signature function defined by 

/:{0,l,...,n-l} — Z 2 = {0,1}, a^f(a) 

and 

| a) = \ai) <g> |a 2 ) ® • • • |a t _i) <g> |a t ), a fc G Z 2 

a = ai2* _1 + a 2 2*~ 2 H h a t _ x 2 + a t , < a < n - 1. 

This operation plays a crucial role in the quantum data-base searching algorithm of 
Grover, |12|] and an important role in quantum computing, in general. Our concern is as 
follows. Is it possible to construct this operator in an efficient manner (steps polynomial 
in t)7 Has such an algorithm been already given in Quantum Computation? 

As far as we know this point is rather unclear, |TT|. See |0 and B. We will discuss 
this point in some detail . 

We would like to construct a road connecting Geometry and Quantum Computation, 
which is not an easy task. We will show one of such attempts as explicitly as possible. 
Though results given in this paper are not new, we do hope our presentation offers new 
perspectives to not only students and/or non-experts but also to experts. 

The contents of this paper are as follows: 

1 Introduction 

2 Grassmann Manifolds 

3 Volume of Unitary Groups 

4 A Question 

5 Quantum Computing 

6 Holonomic Quantum Computation 

Appendix 

A A Family of Flag Manifolds 

B A Generalization of Pauli Matrices 

C General Controlled Unitary Operations 
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2 Grassmann Manifolds 



Let V be a fc-dimensional subspace in C n (0 < k < n). Then it is well-known in Linear 
Algebra that there is only one projection P : C™ — > C n with V = P(C"). Here the 
projection means P 2 = P and Pt = P in M(n; C). 

The Grassmann manifold is in this case defined by all the fc-dimensional subspaces in C n , 
and it is identified with all the projections in M(n; C) with the trace k or the rank k 
(corresponding to V = P(C")). We note that the eigenvalues of a projection are either 
or 1 (by P 2 = P), so the rank of P = trace of P. Therefore we arrive at 

G hn {C) = {P e M(n; C)| P 2 = P, P f = P and trP = k}. (2.1) 

A comment is in order. In general it is not easy to visualize all the fc-dimensional subspaces 
in C n except for experts in Geometry. But it is easy even for us to deal with ( |2.1| ) as will 
be shown in the following. 

We note that G ,n(C) = {0 n } and G nn (C) = {l n }- In particular Gi tTl (C) is called a 
complex projective space and is written as CP n_1 . In (|2.1|) we know a natural symmetry 
(isomorphism) 

k : G k>n (C) — > G n _ fc , n (C), k{P) = l n - P, (2.2) 

so that we have G k>n (C) = G n ^ kjTl {C). 

Now it is easy to see that P can be written as 

P = AE k A~ x for some A e U(n), (2.3) 

where E k is a special projection 

E > - ( o I ) ■ ^ 



n—k 



Therefore we have 
which directly leads to 

In particular 



G k>n (C) = {AE k A~ l \ A e U(n)}, (2.5) 



G ^ s uik^ui-k) • ^ 



G 1>B (C) - CP _ ^^31) " f7(i) " "s 1- ' (2 ' 7) 

see ( ^2|) . Here S fe is the unit sphere in R fc+1 and U(l) = S 1 . We note that G fcjn (C) is a 
complex manifold (moreover, a Kahler manifold) and its complex dimension is kin — k). 

Next let us introduce local coordinates around P in (|2.3|) . We denote by M(n — k, k\ C) 
the set of all in — k) x k - matrices over C and define a map 

V : M{n — k, k; C) — > G k>n (C) 
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as follows 



l fc -zt W l fc o \ ( i k -zt \ 1 



Of course P(0) = P in (p|). 

Here a natural question arises. How many local coordinates do we have on Gfc jn (C) ? The 
number of them is just n C k . 

A comment is in order. We believe that this is the best choice of local coordinates on the 
Grassmann manifold, and this one is called the Oike coordinates in Japan. As far as the 
author knows H. Oike is the first to write down ( |2.8| ), ||19|| . 
From this we can show the curvature form V(Z)dV(Z) A dV(Z): 



V{Z)dV{Z) A dV{Z) 



A 



l fe 

Z ln-fc 



A fc _1 dZt A M n _ k ~ x dZ O 



Z ln-fc 



where 



A fc = l k + Z*Z G M(fc;C) 



n—k 



l n „ k + ZZ ] G M(n-k;C). 



A~\ (2.10) 



(2.11) 



In the following we omit the A symbol and write, for example, VdVdV instead of V(Z)dV(Z)A 
dV{Z) for simplicity. A (global) symplectic 2-form on G kyU {C) is given by 

uj = tiPdVdV 

and its local form 

tr (A- l dZ^M~\dZ) = tr ((l fe + Z^ Z)~ x dZ\\ n - k + ZZ^dZ^j . (2.12) 
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We want to rewrite (|2.12j ). Before doing this let us make some mathematical prelimi- 
naries. For A G M(m, C) and B G M(n, C) a tensor product A<g>B of A and B is defined 

as 



A <g> B = (dijB) for A = (a^) and B = (h 



pqj 



For example, for 



we have 



.4 



0-21 Q>22 



and B 



hi 

&21 i>22 



A®B 



CLl\B dyiB 

a-2\B CL22B 



( Onfeii an&i2 ai2&n a 12 h 12 \ 

an&21 an&22 Oi2&21 &12&22 

021&11 0-21&12 0-22&11 0-22^12 

\ 021^21 a 21^22 a 22^21 a 22^22 / 



(2.13) 
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Therefore componentwise we have (A <g> B)i P j q = AijB pq . Then it is not difficult to see 
tr(A ® B) = tr(A)tr(5), det(A ® B) = {det(A)} n {det(B)} m . (2.14) 

Let us construct a column vector Z in C k ( n ~ k ) from Z in M(n — k, k; C) in a usual 
manner: 



where T means a transpose. Now we rewrite ( [2.12|) as follows: 
W = tr (A, WJICW = tr {d^M~\dZk k l ) 

= i M n~k ® (A^) T } .^dZ^ = (dZ)t { M -_i fc ® (A^f } dZ. (2.15) 

The symplectic volume on G& )n (C) which coincides with the usual volume is given by 

k(n— k) 

{k{n - k)}\ V2 V / ^T 



dv = —, I — ) . (2.16) 



Here 2 ^_ x is a normalization factor. From (|2.15 ) it is easy to see 

Mn-k) = {k{n _ k) y det { M -_l fc g, (A -1)T} Yldz^dZij. (2.17) 



Therefore ( |2.16j ) becomes 



dv = det {M-\ ® (A^f} II (2- 18 ) 



On the other hand by (|2.14j) we have 



det {M~\ ® (A^f } = (detM-_\)* (det^ 1 ))""* = (detM n _ fc )" fe (detA,)"^ 

(2.19) 

Here we note detA^ = detM n _fc. For 



X 

we have 



i fe -zt 

Z l n -fc 



detX = det ( ]t 1 Zt ^ = det ( lfc + Z tZ - = det (l fc + Z f Z) = detA fc . 

V ^ ln-fe / V ^ *-n-k J 



On the other hand 



detX = det ^ ^ ^ J = det ^ ^ x ^ J = det ( ln ~ k + = detM n-^ 



so that 

detAfc = detM n _fc 



From (l2TT9|) det {M~± k ® (A^ 1 ) 7 } = (detA fc ) so we arrive at 



cfo(Z, Zt) = (detA*)"" II ^0 = ( det (^ + Z]Z )Y n II f^=- (2-20) 

From the above-mentioned facts the volume of Grassmann manifold Gr& n (C) is given as 

Vol(G fcn (C))= / 2v^t , 2 21) 

1 , V " y M(n _, )fc;C) {det(U + ZtZ)}" 1 ; 

Problem How can we calculate this integral ? 

3 Volume of Unitary Groups 

Here we will show a heuristic method of evaluation of the volume of unitary group U(n). 
Let S 2k ~ x be the 2k — 1 - dimensional unit sphere (k > 1) over R and the volume be 
Vol(S 2fc_1 ). For example Vo^S 1 ) = 2ir and Vo^S* 3 ) = 2ir 2 . In general we have 

2ir k 

Vol(S^) = (3.1) 



Since we know the fact 



we have 



U(k) 



U(k-l) 



S 2k ~\ (3.2) 



tyW "f/(n-l) X f/(n-2) X X f/(1) X(7(1) 

= S 2 ^ 1 x S 2n " 3 x---xS 3 xS\ (3.3) 

where = means almost equal! 

A comment is in order. Of course the equality does not hold in (|3.3| ) except for the cases 
of n — 1, 2. But for the purpose of volume-counting or coho mo logy-counting there is no 
problem to use ( p.3| )p|. 
By combining ( |3.3| ) and (|3.1|) we obtain 

n n 9 j n(n+l) 

vcp/w) = nvats*-) =n - on.^J-Dr < 3 - 4 > 

Let us evaluate the volume of Grassmann manifold Gfc n (C): 



C/(fc) x C/(n - fc) v K,nv " Vo\(U(k)) x Vol(C/(n - fc)) ' 



2 You will know this "questionable equation" may be rather useful, see for example |2| 
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From 



we obtain 

Vol(G fcift (C)) 



0!1! ••.(* — !)! 0!1! 



[n 



0!1! (n- 1)! 

0!1! ... (A: - 1)! , {n _ k) 



71 



k(n—k) 



(n - k)\ ■ ■ ■ (n - 2)\(n - 1)! 



7T 



(3.5) 



4 A Question 

Combining ( 2.21 ) with ( |3.5|) we have the main result 

0!1! - - - (fc — 1)! 



1-1,3 2y^l 



M(n-k 



k . c) {det(l fc + &Z)} n ( n -k)\---{n- 2)\{n - 1)! 



7T 



k(n—k) 



(4.1) 



It has to be emphasized that the right hand side has been obtained by an indirect path. 
Is it really easy (or practical) to carry out the integration to obtain the right hand side? 
As far as we know, the integral has not been calculated except for the case k — 1. 
Let us review the case k — 1: 



1 



(l + E^k-R™ AA 2 



ndz jdz j 



71 



n-l 



1 



Jj=\ I "J I ) j=l 

The proof is as follows. First let us make a change of variables 
Then we have easily 



(n-l)! 



(4.2) 



(4.3) 



dzjdzj 
2a/^T 



-drAQ 



Under this change of variables ( |4.2| ) becomes 

1 



2tt /"oo 




n_1 ^/fl n_1 



7T 



n-l 



n-l 



'o A (1 + E"=i r i; J= i - i= i 
Here let us once more make a change of variables from (ri, • • • , r ra _i) to • • • , £n-i): 

ri = 6(l"6), 
r 2 = £ik(l-&), 



r n-2 — £l^2 - " " ^n-2(l _ £n-l, 
r n-l = £l^2 " " " ?n-2Sn-l- 
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Conversely we have 



6 = n +r 2 H hr n _2 + r n _i, 

r 2 H h + r n _i 



n-2 



T\ + r 2 H h r n _ 2 + r„_i 

7"n-2 + rV-1 
+ ^n-2 + ?"n-l ' 



Sn— 1 ) 

r n -2 + r„_i 
0<£i<oo, 0<6,-"-,a-i < 1, 

n—1 n—1 

n^=^ n "V~ 3 ---a-2n^- 

3=1 3=1 

Under this change of variables (|4.2|) becomes 



/oo c n—2 />1 />1 

(l + 6 )n ^ / ^ n " 3 ^2 • • • J £n-2^„-2 

-1 /-l /-l 

n-1 / £ n— 2 j<- I (■ n— 3, 



= 7r «-l, 6 n ~^6---/ fn-2#„-2 

Jo Jo JO 

= n n - 1 — - = 7 ^ ■ (4.4) 

n- In- 2 2 (n-1)! v ; 

The direct proof of ( |4.1|) for k = 1 is relatively easy as shown above. But for > 2 we do 
not know such a proof (a direct proof may be very complicated). Therefore let us present 



Problem Give a direct proof to 

/ 



M(n-*, fc; c) {det(l, + ZtZ)}" (n - A;)! • • • (n - 2)!(n - 1)! 



As for another approach to the above problem we refer to JRJ. In this paper coherent 
states based on Grassmann manifolds have been constructed. 



5 Quantum Computing 

Let us move to the main subject of Quantum Computing. The typical examples of 
quantum algorithms up to now are 



Factoring algorithm of integers by P. Shor f28| 



Quantum data-base searching algorithm by L. Grover 
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See and P5| , or [nj for general introduction. [TJ[] and |T5| are also recommended. 

In Quantum Computing we in general expect an exponential speedup compared to 
classical ones, so we must construct necessarily unitary matrices in U(n) in an efficient 
manner when n is a huge number like 2 100 . 

Problem How can we construct unitary matrices in an efficient manner? 

Let us come back to fl2~T|). We denote the set ofnxn projection operators by 

G n (C) = {Pe M(n; C)| P 2 = P, P ] = P}. (5.1) 

The elements of G n (C) are classified by the trace, so G n (C) can be decomposed into a 
disjoint union 

n 

G n (C)= \jG k>n (C). (5.2) 

k=0 

For a fc-dimensional subspace V in C n (0 < k < n) let {vi, V2, • • • , v^} be an orthonormal 
basis (namely, < Vj, v,,- >= <5y) and set 

V r = (v 1 ,v 2 ,---,v fc )eM(n,fc;C). (5.3) 

We have identified a /c-dimensional subspace V with a matrix V in ( |5.3| ) for simplicity 
(maybe there is no confusion). Then we have an equivalence 

{ v i> v 2> ••• ; v fc} : orthonormal V'V = lfc. 

Then it is easy to see that all orthonormal basis in are given by 

{Va\ a G U(k)}. (5.4) 

The projection corresponding to V is written by 

p = yyt eG M (c). (5.5) 

We remark that (Va)(Va)T = VaasV' = VV* = P, namely P is of course independent of 
a G U(k). This P is also expressed as 

k 

P = E V ^V- (5-6) 
If we use Dirac's bra-ket notation Vj = \j), then P = Yuj=i ^his notation may be 



popular in Physics rather than (]5.6|) . 

How can we construct an element of unitary group from an element of Grassmann 
manifolds? We have a canonical method, namely 

G n (C) — > C/(n) : P ^ C/ = 1„ - 2P (5.7) 
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This U is called a uniton in the field of harmonic maps. Moreover we can consider a 
product of some unitons, namely, for any S C {0, 1, • • • , n — 1, n} 



U = H{ln - 2P S ) for P j e G jtn (C) 



(5. 



In particular 



ra-1 



U = \[{l n -2P j ) for P j e G iin (C) (5.9) 

is very important in the field of harmonic maps, see for example |4[] and J3TJ . 
Many important unitary matrices are made by this wayQ. 

These unitary matrices also play an important role in Quantum Computing as shown in 
the following. 

We consider a qubit (quantum bit) space of quantum particles. The 1-qubit space is 
identified with C 2 with basis {|0), |1)} ; 



Vect c {|0),|l)}, |0> 



|1> 



The qubit space of t-particles is the tensor product (not direct sum!) of C 2 



with basis 

For example, 
|0>®|0> = 



{|ni,n 2 , ...,n t ) = |m) ® \n 2 ) 



C 2 = (C 



2\®t 



(5.10) 



\n t ) | rij ez 2 = {0,1} } 



/ 1 \ 




v o y 



, |0)®|1) 



1 



v o y 



|1) ® |0) 





1 

v o y 



, ll>® II) 






V 1 / 



Now we take the Walsh- Hadamard transformation W defined by 

W:\0)-* -L(|0> + |1», W:|1> — -L(|0>-|1», (5.11) 



in matrix notation, 



W=±(\ _\ \ G 0(2) c 17(2). (5.12) 



3 Those used in j2| for data-base searching algorithms are of this form with appropriate Pj 
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This transformation (or matrix) is unitary and it plays a very important role in Quantum 
Computing. Moreover is easy to realize it in Quantum Optics. Let us list some important 
properties of W: 



W 2 = 1 2 , W ] = W = W~ 
a x = Wa 3 W~\ 



where {<7i, <t 2 , a 3 } are the Pauli matrices: 

1 



1 



(5.13) 
(5.14) 



(5.15) 



See Appendix B for a generalization of Pauli matrices. Next we consider t-tensor product 
of W (t e N) : 

W® 1 = W <g> W <g> • • • ® W (t-times) . (5.16) 
This matrix of course operates on the space (|5.10 ). For example 



and 



W <g> W ® W = —= 





/ 1 


1 


1 


1 ^ 








1 


1 


-1 


1 


-1 








~ 2 


1 


1 


-1 


-1 










V 1 


-1 


-1 


1 J 








/ 1 


1 


1 


1 


1 


1 


1 


1 \ 


1 - 


-1 


1 


-1 


1 - 


1 


1 


-1 


1 


1 


-1 


-1 


1 


1 


-1 


-1 


1 - 


-1 


-1 


1 


1 - 


1 


-1 


1 


1 


1 


1 


1 


-1 - 


1 


-1 


-1 


1 - 


-1 


1 


-1 


-1 


1 


-1 


1 


1 


1 


-1 


-1 


-1 - 


1 


1 


1 


u - 


-1 


-1 


1 


-1 


1 


1 


-1 J 



(5.17) 



(5.18) 



Hereafter we set n = 2*. Then ( |5.16|) means W® 1 G U(n). The very important fact is 
that ( |5.16| ) can be constructed by only t{= log 2 (n))-steps in Quantum Computing. Let 
us show the matrix-component of (|5.16|) is given by 

(h, i 2 , . . .,i t \W & \j u j2, ...Jt) = -^(-l)£Ui*i* (5.19) 

or if we set 

|i) = |ii)®|ia)®---®|it), i = i^ 1 ' 1 + i 2 2^ 2 + ■ ■ ■ + i t , 0<i<n-l 
we have 

(i\W®t\j) = 4=(-l)*" i (5-20) 

in 



where % ■ j means the sum of bit- wise products Ylk=i ikjk- 
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The proof goes as follows. From ( |5.12| ) we know 

w\i) = -J=(io) + (-i) 4 |i» = -L E(-ir^) ; 

which implies 

W^\j) = W\h) ® W\j 2 ) ® ■ • • © 

4 E E ■ • • ^(-i)^ 1 ^^-^!^) (8 1*5,)®... |fc t ) 



fc 1= 0fc 2 =0 k t =0 
n-l 



j^J2(-ir\k). (5.21) 



fc=0 

Therefore we obtain 

n— 1 -, n—1 



1 " 1 1 

= 7^ E(- 1 ) fe "(^) = 7^ E(- 1 )^^ = " 

V fe=0 ^ fe=0 v 

Moreover ( 5.16 ) has an interesting property which we can guess from (|5.17 ) and ( 5.18| ): 



n-l 



3=0 



or 

n-l 



= { ^ others. < 5 ' 23 ) 

Let us clarify the meaning of ( 5.20|) from the point of view of Group Theory. 
We note that Z 2 is an abelian group with operation © 

© = 0, © 1 = 1, 1 e = 1, 1 © 1 = 0. (5.24) 
Then Z 2 * is a natural product group of Z 2 . We denote its element by 
i = (zi, zi, • • • , i t ) < — ► i = i x 2 l ~ l + i 2 2*~ 2 + --- + i t . 
For i e Z 2 * we define 

Xi : Z 2 * — > C* = C - {0}, Xi(j) = v^l^lj) = (-1)^'. (5.25) 
Then we can show that 

Xi(j©k) = Xi (j) Xi (k). (5.26) 

That is, Xi is a character of the abelian group Z 2 *. 
The proof is as follows. From ( |5.24j ) we know 

x © y = x + y — 2xy for x,yGZ 2 . (5.27) 
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Therefore we obtain 



Xi (j©k) = (— l)^*=i n(n® k i) = (_i)E?=i i iO'i+fe i -2i i fc i ) _ r_i\Ei=ih{h+ki) 

= (_i)Ef= a ni i (_ 1 )Ef= 1 ^ ! = = Xi (j) Xi (k) ■ (5.28) 



These characters play an important role in Discrete Fourier Transform, see |14] or recent 
preprint [IS . 

Now we consider a controlled NOT operation (gate) which we will denote by C-NOT 
in the following. It is defined by 



(5.29) 
(5.30) 



C-NOT : 


|0,0> 


-|o,o), 


|o,i) 


- |o,i) 




|i, o> 


->|M), 


|M> 


^ |i,o) 


or more compactly 










C-NOT : 


a, b) — 


— > a,o(E 


6), a 




Graphically it is expressed as 










a) 




-• 




a) 



.X a \b) = \a@b) 



and the matrix representation is 



C-NOT 



/ 1 \ 

10 

1 

\ 1 J 



(5.31) 



Here a@b = a + b (mod 2) and we note the relation 

X a \b) = al\b) = \a®b) for a, b G Z 2 . 



(5.32) 



In this case the first bit is called a control bit and the second a target bit. 

Of course we can consider the reverse case. Namely, the first bit is a target one and the 

second a control one, which is also called the controlled NOT operation: 



C-NOT : \a,b) 
and the matrix representation is 

C-NOT 



\a © b, b), a, b E Z 2 , 



/ 1 \ 

1 

10 

\ 1 / 



(5.33) 



(5.34) 
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A comment is in order. In the 1-qubit case we may assume that we can construct all 
unitary operations in U(2) (we call the operation universal). In the 2-qubit case how 



can we construct all unitary operations in U{4)7 If we can construct the C-NOT ( 5.31 ) 
( |5.34| ) in our system, then we can show the operation is universal, see [[| and [[]J . This is 



a crucial point in quantum computing. We comment here that the C-NOT ( 5.31 ) can be 
written clS cL uniton ( p.7| ) 

(5.35) 



where 



C-NOT 




1 4 - 


2P 






(0 











\ 


1 
















2 








1 


-1 











-1 


1 


/ 



(5.36) 



and this P can be diagonalized by making use of Walsh- Hadamard transformation ( 5.12|) 
like 



P = (1 2 ® W)E 1 (1 2 <g> Wy 1 = (1 2 <g> W)(a 1 ® a 1 )E 1 {a l ® «7i) _1 (l2 ® W^ 1 (5.37) 
where 



/0 



V 











1 / 



( 1 



V 











0/ 



(5.38) 



More generally for the t-qubit case we can construct (t — l)-repeated controlled- not 
operator and show it is a uniton. 

The (t — l)-repeated controlled-not operation is defined by 



C (t - 1} -NOT : lax, 



a>2, 



<k-i,<h) 



ai, 02, 



, a t _i, aia 2 ■ • -a t _i © a*), 



or in matrix form 



/ 1 



C^-NOT 



a k G Z 2 (A; = 1,2, •••,*), 



(5.39) 



V 




1 



1 

0/ 



2* x 2* - matrix. 



(5.40) 



As for the explicit construction of (t — l)-repeated controlled-not operator see [|TJ or 
0. See also Appendix. But unfortunately the construction is not efficient ! 
A comment is in order. In |l] a rough estimation of the number of steps to construct the 
operator ( f5.39| ) is given and it is confirmed that an efficient construction is possible. But 
no explicit construction is given. 
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By the way, since 



/ W 



1 2 ®(*-i) iy 



W 



we have 



W / 
/ 1 



(5.41) 



(la®^- 1 ) ^ w) C^-NOT (l 5 



9(t-l) 



w) 



\ 



1 



V 



1 / 



= l n -2|n-l)(n-l|. (5.42) 

Therefore the construction of C^^-NOT and l n — 2\n — l)(n — 1| have almost the same 

number of steps. Is it possible to construct this operator efficiently ? 

As far as we know an explicit and efficient construction of this operation has not yet been 

givenf]. 



Problem Give an explicit and efficient algorithm to this operation. 
Let us consider a set 



where 



Fk — In. — 2E, 



{F k | 1 < k < n} 



(5.43) 



<-n—k 



If Fx can be constructed, then the other F's can be easily obtained. First let us show this 
with a simple example (t=2) : 



Fx 



( ~ l 
\ 



\ 



1 / 



1 4 -2|0)(0|. 



Now we set 



Ux = l 2 ® ox 



U 3 = ox® ox 



( 1 
1 

V 



V 1 



\ 

1 

1 \ 



/ 



U 2 = ox® 1 2 



1 \ 

1 



l I have not yet succeeded in such a construction. 
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then we have 



U 1 F 1 U 1 



U 2 F X U 2 



U 3 F,U 3 = 



( 1 



/ 1 



/ 1 



V 



1 / 



- 1 / 



l4-2|l)(l|, 



1 4 -2|2)(2|, 



= 1 4 -2|3)(3|, 



so that it is easy to check 

F 1 (U 1 F 1 U 1 ) = F 2 , F 2 (U 2 F 1 U 2 ) = F 3 and F 3 (U 3 F 1 U 3 ) = -1 4 . (5.44) 
Let us prove the general case. For i = + i 2 2 f ~ 2 + • ■ ■ + i t (0 < i < n — 1) we set 



Since 



Ui = al 1 <g> a? <g> • • • <g> <r«, (Uj = U l = Ur 1 ). 



Ui\0) = al 1 <g> al 2 ■ ■ ■ crj* (|0) <g> |0) <g> • • • <g> |0)) 
= <t\ l \0)®<t?\0)®---®<t?\0) 
= \ii) <g> \i 2 ) <g> ■ ■ ■ <g> \U) = 



(5.45) 



we have 



l n - 2\i)(i\ = ^ (l n - 2|0)(0|) Ui = UiF.Ui 
Therefore it is easy to see that 

F k {U k F l U k ) = F k+l {\<k<n-l). ■ 
We note that this procedure is not efficient. 



(5.46) 
(5.47) 

(5.48) 



Now, let us make a comment on Grover's data-base searching algorithm. In his algo- 
rithm the following two unitary operations play an essential role: 



l n -2|i)(i|, l n -2|s)(s|, 
in which the state \s) (s stands for sum) is defined by 



n— 1 

*) = -^£lO = w®*|o). 



(5.49) 



(5.50) 
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We find via (|5"2IQ that 

l n - 2\s)(s\ = W m (l n - 2|0)(0|) W m = W m FxW m . 



(5.51) 



Namely, the two operations ( 5.49Q are both unitons and can be diagonalized by the efficient 
unitary operations Ui and W® 1 . 

Finally, let us mention about a relation between (t — l)-repeated controlled- not op- 
eration and F\. Since 



2\n- l)(n- 1| = Un^FiU, 



n—l 



(l^*- 1 ) ® W) C ( ' _1) -NOT (1 2 ®(*- 1 ) <g> W) = 1 

by fl5.42| ), we have 

Fx = U n -! (1 2 0( *- 1} ® W) C ( '~ 1} -NOT (la®^ 1 ) ® W) C/ n _i. 

By substituting C/„_i = oi (g> cri ® • • • <S> o"i ® cri into the equation above we arrive at the 
desired relation 



Ft = (ai 



O(t-l) 



triW) C^-NOT (a^^ 1 ) ® Wo-i) 



(5.52) 



6 Holonomic Quantum Computation 

In this section we briefly introduce a simplified version of Holonomic Quantum Computa- 
tion. The full story would require detailed knowledge of Quantum Mechanics, Quantum 
Optics and Global Analysis, so it will have to wait for another occasion. 



This model was proposed by Zanardi and Rasetti |3T[ and |22] and it has been devel- 
oped by Fujii §, §, 0, § and Pachos 0], |21|. 

This model uses the non-abelian Berry phase (quantum holonomy in the mathematical 
terminology [0) in the process of quantum computing. In this model a Hamiltonian (in- 
cluding some parameters) must have certain degeneracy because an adiabatic connection 
(the non-abelian Berry connection) is introduced in terms of the degeneracy, see [27|. In 
other words, a quantum computational bundle is introduced on some parameter space 
due to this degeneracy and the canonical connection of this bundle is just the one above. 
On this bundle Holonomic Quantum Computation is performed by making use of the 
holonomy operations. We note our method is completely geometrical. 



Here we introduce quantum computational bundles, ||, 0, and ||. For this pur- 
pose we need universal principal and vector bundles over infinite dimensional Grassmann 
manifolds. We also need an infinite dimensional vector space called a Hilbert (or Fock) 
space. 

Let H be a separable Hilbert space over C. For m G N, we set 

st m {H) = {v = ( Vl , ■ ■ ■ , Vm ) e n x ■ • • x n\ vW = i m } , (6.i) 

where l m is a unit matrix in M(m, C). This is called a (universal) Stiefel manifold. Note 
that the unitary group U(m) acts on St m {7i) from the right: 

St m (H) x U{m) -> St m (H) : (V, a) i-> Va. (6.2) 
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Next we define a (universal) Grassmann manifold 

Gr m (H) = {Xe M(H)\ X 2 = X,X ] = X and trX = m) , 



(6.3) 



where M(TC) denotes a space of all bounded linear operators on 7i. Then we have a 
projection 

m 

vr : St m (H) - Gr m (H), tt(V) = VV^ = ^v jV }, (6.4) 

j=i 

compatible with the action Q (ir{Va) = Va{Va) ] = Vaa^V^ = VV ] = 7r(V)). 
Now the set 

{f/(m),5t m (H),7r,Gr m (H)}, 
is called a (universal) principal ?7(m) bundle, see fl8j and @. We set 



E m {H) = {(X,v) e Gr m {H) xH\Xv = v}. 
Then we have also a projection 

vr : E m (H) -> Gr m (H) , tt((X, «)) = X 

The set 



(6.5) 
(6.6) 



(6.7) 
(6.8) 



{C m ,E m (H),7r,Gr m (H)} 

is called a (universal) m-th vector bundle. This vector bundle is associated with the 
principal U(m) bundle (|6.5| ). 

Next let M be a finite or infinite dimensional differentiable manifold and the map 
P : M Gr m (7i) be given (called a projector). Using this P we can define the pullback 
bundles over M from ( |6.5|) and (|6.8j ): 

U(m),St, n Jv M} = P* {U(m), St m (H), vr, Gr m (H)} , (6.9) 

{c m ,E,7r^,M} =P*{C m ,E m (n),7r,Gr m (n)}, (6.10) 

see [[ill. Of course the second bundle ( |6.10|) is a vector bundle associated with the first 
one (p|): 

U(m) U(m) 



St 



M 



St m {H) 



Gr m (H) 



C 



E 



M 



E„ 



(H) 



Gr m {H) 



Let M. be a parameter space (a complex manifold in general) and we denote by A its 
element. Let Ao be a fixed reference point of M.. Let H\ be a family of Hamiltonians 
parameterized by M. acting on the Fock space Ti. We set H = H\ for simplicity and 
assume that this has an m-fold degenerate vacuum: 



H Vj = 0, j = l,...,m. 



(6.11) 
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These v^s form an m-dimensional vector space. We may assume that (vi\vj) = 5ij. Then 
Oi, • • • ,v m ) e St m (H) and 

Namely, F is a vector space associated with the o.n. basis (vi, ■ • ■ , v m ). 

Next we assume for simplicity that a family of unitary operators parameterized by A4 

W U(H), W(X ) = identity, (6.12) 

connects H\ and H isospectrally: 

H x = W(X)H W(X)- 1 . (6.13) 

In this case there is no level crossing of eigenvalues. Making use of W(X) we can define a 
projector 

P : M - Gr m (H), P(A) = W{X) [fl v 3 v jj W W' 1 ( 6 - 14 ) 
and the pullback bundles over Ai: 

{u(m),St, tt Jv m} , {C m , E, tt^m}. (6.15) 

For the latter we set 

\vac) = (vi, ■■■,v m ). (6.16) 

In this case a canonical connection form A of the principal bundle j C/ (m) , St, 7Tg^, .M j is 
given by 

^4 = ( wac |iy(A)- 1 c /iy(A)|t;ac), (6.17) 
where d is a differential form on M. 



together with its curvature form (see [27 and [IS]) 

F = dA + A/\A. (6.18) 

Let 7 be a loop in M. at A , 7 : [0, 1] — > M, 7(0) = 7(1) = A . For this 7 a holonomy 
operator T_4 is defined: 

1^(7) = Vexp |y .4 j G I7(m), (6.19) 

where P means path-ordering, see for example |H] . This acts on the fiber Fo at Ao of the 
vector bundle |c m , E, n^, m\ as follows: x — ► T_4(7)x. The holonomy group Hoi (A) 
is in general a subgroup of U{m) . In the case of Hoi {A) = U(m), A is called irreducible. 
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The irreducibility of A is very important because it means the universality of quantum 
computation. To check whether A is irreducible or not we need its curvature form ( |6.18|) , 

see 



In the Holonomic Quantum Computation we take 



Encoding of Information = 
Processing of Information 



x€ F , 

> r^(7) : x -> T^(7)x = Ax. 



(6.20) 



See the following figure. 




Our model is relatively complicated compared to the other geometric models and much 
more so than the usual spin models. We have a lot of problems to solve in the near future. 
We strongly hope young graduate students to invigorate this field. 

Acknowledgment. The author wishes to thank Yoshinori Machida for his warm hospitality 
at Numazu College of Technology. I also wishes to thank Akira Asada, Ryu Sasaki and 
Tatsuo Suzuki for reading this manuscript and making some useful comments. 



Appendix 



A. A Family of Flag Manifolds 

Let us make a comment on an interesting relation between flag manifolds and the kernel 
of the exponential map defined on matrices. Here a (generalized) flag manifold (which is 
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a useful manifold as shown in the following) is a natural generalization of the Grassmann 
one. 

First of all we make a brief review. For 



2-nJ-lt 



exp : R — > S 1 C C, exp(t) = e 

the kernel of this map is ker(exp)=Z C R. 

We define by H(n, C) the set of all hermitian matrices 

H(n, C) = {X G M(n, C) | X f = X}. 

Of course H(l, C) = R. Note that each element of H{n, C) can be diagonalized by some 
unitary matrix. 

The exponential map is now defined as 

E : H{n, C) — ► U(n), E(X) = e 27TV ^ x . (A.l) 
Here our target is ker(E). 
Problem What is the structure of ker(E) ? 

Our claim is that ker(E) is a family of flag manifolds. For that we write ker(E) as 

K n (C) = {Xe H(n, C) | e 2 ^ ix = l n }. (A.2) 

First we prove 

G n (C) C K n (C). (A.3) 
Because since P 2 = P from the definition, P k = P for k > 1, so that 

k=l ' k=l 

(A.4) 

We will prove that G n (C) becomes a kind of basis for K n (C). 

For X G K n (C) we write the set of all eigenvalues of X as spec(X). Then spec(X) = 
{0,l}forXeG„(C). 

It is clear that spec(A) C Z. For X e K n (C) we have 

3 

spec(X) — {rii(di), rik(d k ), ■■■,rij(dj)} where G Z and ^^rffc = n, (A. 5) 

k=i 

in which (dfe) is the multiplicity of the eigenvalue n^. Since X is diagonalized by some 
U G U(n), 

j 

X = UX U- 1 = Y^n k P dk , (A.6) 
k=i 
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where 



( n{± ix 



\ 



n 2 l 



\ 



( 



di 



\ 



u 



1<& 



V 



(A.7) 



Here we list some properties of the set of projections {P dk }'- 

(1) P dk e G dfe , n (C), (2) P dk P dl = 5 kl P dl , (3) P dl + P d2 + ■ ■ ■ + P d] = l n . 

Let us here prepare a terminology. For X G -K'n(C) we call set of the eigenvalues together 
with multiplicities 

{(ni, di), (n 2 , d 2 ), • • • , {rij, dj)} 

the spectral type of X. 

Then it is easy to see that X and Y G fT n (C) are of the same spectral type (X ~ Y) if 
and only if Y = UXU~ l for some U G U(n). For X G X n (C) we define 

c(x) = {y g # n (c) | y ~ x} . 

We have clearly C(X) = C(Xq). Then it is easy to see that K n (C) can be classified by 
the spectral type 

K n (C) = |J C(X) = |J C7(X ) (A.8) 

and the unitary group Z7(n) acts on C(X) as follows: 

U(n)xC(X) — >C(X) : (U, X) ^ UXU' 1 . 
Since this action is free and transitive, the isotropy group at X is 

U{di) x U(d 2 ) x ■ ■ • x C/(^) , 

so that we have 

U(n) 



C(X) 



U{d x ) x U(d 2 ) x • • • x ^(dj) ' 

The right hand side is called a generalized flag manifold. In particular when d\ = d 2 
■ • ■ — d n — 1 (there is no overlapping in the eigenvalues of X) we have 

U(n) 



(A.9) 



C(X) 



(A.10) 



U(l) x U(l) x ••• x U(l) ' 
This is called a flag manifold. 

Namely by ( |A.8|) we know that K n (C) is a family of generalized flag manifolds. 

A comment is in order. For the Grassmann manifolds we have very good local coordi- 
nates like (|2.8|) , while we don't know good local coordinates for generalized flag manifolds. 

Problem Find a good local coordinate system. 



For some applications of generalized flag manifolds the paper |23| is recommended. See 
also |^3| and references therein. 
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B. A Generalization of Pauli Matrices 



Here let us introduce a generalization of Pauli matrices ( |5.15|) which has been used in 
several situations in both Quantum Field Theory and Quantum Computation. 

First of all we summarize the properties of Pauli matrices. By fl5.15| ) cr 2 = V — lo"iO"3, 
so that the essential elements of Pauli matrices are {cxi, 03} and they satisfy 



cr; 



(To 



L2- 



t 



O x = CTi, CT 3 = (J3 



a 3 a 1 



-0-1(73. 



(B.l) 



Let {Si, S3} be the following matrices in M(n, C) 

1 \ /I 



Si 



/ 

1 
1 

1 



V 



1 oj 



a 



where a is a primitive root of unity a 1 

n-l 



a 



1 ( cr = e * ). We note that 



n-l 



(B.2) 



a = a 



n-l 



0. 



The two matrices S 3 } are generalizations of Pauli matrices {01, 03}, but they are not 
hermitian. Here we list some of their important properties: 



V-171— 1 

^3 



En yin -| vm y-in— 1 y-if 

1 _ Zj 3 — L n , Zj 1 — Zj 1 , Zj 3 

If we define a Vandermonde matrix W based on a as 



S3S1 — (JS1S3 



(B.3) 



/ 1 1 

1 a 11 ' 1 
1 o n ~ 2 



W 



1 

a 2(n-l) 



a 



1 

(n-l) 2 



cr 



2(n-2) 



cr 



(n-l)(n-2) 



1 

V 1 

/ 1 
1 
1 



w* = A= 



cr" 
a 

1 

cr 
a 2 



a 2 



CT 



a 



2(n-l) 



(7 



n-l 
1 

a-n-i 

a 2(n-l) 



(B.4) 



1 Cr™" 2 
,n-l 



cr 



V 1 ^ 



cr 



2(n-2) 
2 (n-l) 



a 



(n-l)(n-2) 



cr 



(B.5) 



then it is not difficult to see 



£1 = WE 3 W^ = WE 3 W-\ 



:b.6) 
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For example, for n = 3 




Si, 



where we have used that cr 3 = 1, a = a 2 and 1 + a + a 2 = 0. 
That is, Ei can be diagonalized by making use of W. 

A comment is in order. Since W corresponds to the Walsh- Hadamard matrix (|5.12 ), so 
it may be possible to call W the generalized Walsh- Hadamard matrix. 



C. General Controlled Unitary Operations 

Here let us introduce a usual construction of general controlled unitary operations to help 
the understanding of general controlled NOT one. In the following arguments if we take 
U = X = o\ then they reduce to the arguments of a construction of general controlled 
NOT operator. 

First of all let us recall ( |5.27| ). For x,y,z G Z 2 we have identities : 

x + y — x © y = 2xy, (C.l) 
x + y + z — x ®y — x ® z — y®z + x®y®z = 4xyz, (C2) 

where x ® y = x + y (mod 2). For the most general identities of above type see U and 



The controlled-controlled unitary operations are constructed by making use of both several 
controlled unitary operations and controlled NOT operations: Let U be an arbitrarily 
unitary matrix in £7(2) and V a unitary one in U(2) satisfying V 2 = U . Then by ( |C.1| ) 
we have 



yx+y-x®y _ ylxy _ ^y2\xy _ jjxy 

yx+y-x(By _ yxyyy-x®y _ yxyy ry-l\x®y _ yxyy ryUx®y 

so a controlled-controlled U operation is graphically represented as 



(C.3) 
(C.4) 




U xy \z) 
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You should read this figure as follows: From the left to the right 

\x) <g> \y) © \z) — ► \x) © \y) © U xy \z). 

The controlled-controlled-controlled unitary operations are constructed as follows: Let 
U be an arbitrarily unitary matrix in U(2) and V be a unitary one in £7(2) satisfying 
V 4 = U. Then by Q 

^i+i/+2-(a:0i/+3;92+i(ffiz)+3;ffi!/ffi2 _ yixyx _ jjxyz 

yrr+?/+z-(iEffij/+:rffiz+j/ffiz)+:Effij/ffiz _ yxyyyz ny\\x®y /y\\X®z /yf \3/ffizyc©2/©z /q g\ 

so a controlled-controlled-controlled C/ operation is graphically represented as 



XHX 



x x 



30 l*> 



k) — ®-(£h2) — ® UTU » 

This figure means 

\x) © \y) © \z) © \w) — ► \x) © \y) © \z) © U xyz \w). 
A comment is in order. For the case U = X = a\ we have from ( p. 39 ) 

|<2i, 02, • • • , a n -l, a>n) H a l> a 2i • • • , fln-lj ' ' ' Qn-1 © O n ) 

=|oi) © I02) © ■ • • |a n _i) © |oia2 ■ • • a n _i © a„) 

=|ai> © |a 2 ) © • • ■ |a n _x> © X s1 "'"*- 1 |a n ). (C.7) 

As can be seen from the figures the well-known construction of general controlled 
unitary operations needs exponential steps. Namely it is not efficient. For more details 
see and 0. 
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